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00 . Abstract 

Q ■ We study the A^j-flavor Gross- Neveu model in 2+1 dimensions with a baryon chemical 

potential /i, using both analytical and numerical methods. In particular, we study the 

^ ' self-consistent Boltzmann equation in the Fermi liquid framework using the quasipar- 

ticle interaction calculated to 0(1/ Nf), and find solutions for zero sound propagation 
for almost all /i > fie, the critical chemical potential for chiral symmetry restoration. 
O-i, Next we present results of a numerical lattice simulation, examining temporal cor- 

^ ■ relation functions of mesons defined using a point-split interpolating operator, and 

finding evidence for phonon-like behaviour characterised by a linear dispersion rela- 

/\t ' tion in the long wavelength limit. We argue that our results provide the first evidence 

c^ ! for a collective excitation in a lattice simulation. 

PACS: ll.lO.Kk, ll.15.Ha, ll.15.Pg, 21.65.+f, 67.80Cx 
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1 Introduction 

Such is the perceived difficuhy of performing non-perturbative lattice simulations of 
systems with no n- zero baryon chemical potential /i that it is often overlooked that, 
in those cases where simulations are possible, the information extracted considerably 
exceeds that coming from better understood simulations at non-zero temperature 
T. The reason is that for ^/T ^ 1 there is no stringent limitation on the number 
of temporal lattice spacings Lj, so that correlation functions in Euclidean time can 
be studied with good resolution. These functions encode all information about the 
spectrum of the system under study, including energies and widths of all quasiparticle 
and collective excitations. For instance, in Two Color QCD spectral studies have 
revealed a scalar isoscalar Goldstone boson signalling the superfluid nature of the high 
density state IJ, and evidence for an in- medium decrease in mass of the p- meson j2j. 
In the NJL model in 2+ld the quasiparticle dispersion relation E{k) is smooth for 
k ':^ kp implying a well-defined Fermi surface, but anomalously slow temporal decay 
in the scalar isoscalar diquark channel signals a critical state characteristic of a thin 
film superfluid jHlll]. In NJLs+i by contrast, the scalar isoscalar diquark channel has 
a Goldstone pole, and the quasiparticle dispersion is gapped, providing evidence for 
superfluidity via orthodox BCS diquark pair condensation [H]. 

In refererence jH] we have studied a still simpler model, the 2+ld Gross-Neveu 
(GN) model describing Nf flavors of self- interacting fermions, with continuum La- 
grangian 

C = i,[(^ + ^^,)ij - ^{i,^f . (1) 

As /i is increased, the model exhibits a sharp flrst-order transition at a critical ^c 
from a phase where a fermion mass Sq is dynamically generated but baryon density 
fT-B = {4'loi') vanishes to one where Sq = and n^ oc /i^ [7j. By examining various 
temporal correlation functions we were able to demonstrate many features of this 
high density phase. For example, the quasiparticle dispersion relation is again con- 
sistent with a sharp Fermi surface, a feature reinforced by the oscillatory behaviour 
of meson wavefunctions as the spatial separation between quark and antiquark flelds 
is increased, reminiscent of a phenomenon in many body physics known as Friedel 
oscillation. Secondly, we studied the correlator of the scalar meson a and found that, 
in contrast to its behaviour in vacuum, in the high density phase it decays with a pole 
mass Mo- ~ 2\J fi^fi — fj^c), the value predicted to leading non-trivial order in 1/Nf. 
Finally, we studied meson correlators at non-zero momentum with several different 



quantum numbers, and found that in general they decay algebraically rather than 
exponentially with Euclidean time, since it is usually possible to excite a particle - 
hole pair at the Fermi surface with zero energy cost. 

There is a phenomenological model of interacting degenerate systems, originally 
due to Landau jH] , known as the Fermi liquid. From a particle physics perspective it is 
conceptually simple; the dominant low energy excitations, the so-called quasiparticles, 
carry the same quantum numbers (spin, isospin, baryon number, . . . ) as free particles 
and holes, and can be viewed simply as dressed versions of the elementary quanta. 
Beyond the explicit breaking of Lorentz symmetry introduced hj fi ^ ^ ub y^ 0, 
there is no broken symmetry or order parameter which distinguishes the Fermi liquid 
from non-interacting degenerate fermions. However if certain conditions are fulfilled, 
such as an effective interaction between quasiparticles which is both short-ranged 
and repulsive, then there is a massless bosonic excitation in the spectrum as T ^ 0. 
Carrying zero baryon charge, it is a phonon, ie. a quantum of a collective excitation 
called zero sound [9 . Sound propagation occurs in any elastic medium; zero sound 
happens when the elasticity originates not from collisions between individual particles, 
but from the force on a single particle due to its coherent interaction with all others 
present in the medium. It can be pictured as a propagating distortion in the local 
shape of the Fermi surface, whose speed exceeds that of conventional "first" sound. 

In jH] we found one meson channel whose temporal decay resembled that of an 
isolated pole rather than that of a particle - hole continuum, and which moreover 
yielded a phonon-like dispersion relation u oc k. In the current paper we seek to 
build on our tentative identification of this state with zero sound, firstly by finding 
an analytic solution appropriate for GN2+1, and secondly by performing a more refined 
numerical analysis in the channel of interest. In Sec. |21we review the derivation of 
the self-consistent equation for zero sound propagation in the framework of orthodox 
Fermi liquid theory, and in Sec. El we solve this equation for the particular case of 
GN2+1, using the 0{1/Nf) expression for quasiparticle interactions derived in |6.. We 
find solutions for sound propagation for almost all /i > fie, with speed Po > Pf, the 
Fermi velocity, which is the characteristic scale for a degenerate system. In Sec. 0] 
we present results of a numerical simulation of GN2+1; we apply orthodox meson 
spectroscopy techniques and find evidence for a state with dispersion u{k) oc k in the 
small-fc limit. The propagation speed is of the same order, but slightly less than, the 
Fermi velocity. Our results are summarised with a discussion about possible sources 
for the discrepancy in Sec. El 



2 Theoretical Background 

We begin with a brief introduction of the theory behind sound propagation in de- 
generate systems, leaning heavily on the treatment in JU]- Consider a system of 
(quasi)particles with phase space distribution n given by 

n{k,x,t) = nQ{k) + 6n{k,x,t) (2) 

and energy 

e{k,r,t) = eo{k) + 6e{k,f,t), (3) 

where for states near the Fermi surface the equilibrium conditions are given by 

no(^) = (exp(£o(^)-/i)/T+l)-^ 

6oik) ~ fi + (3pi\k\-kp). (4) 

The parameters are Fermi momentum kp, and Fermi velocity (3f = \^k£o{\k\ = A^f)!, 
which in general differ from their free field theory values. Small departures from 
equilibrium are described by 

Se{k,x,t) =tT / J^{k,k')Sn{k',x,t), (5) 

where d is the number of spatial dimensions and the trace is over internal degrees of 
freedom such as spin and flavor. The Fermi liquid interaction J-" encodes the response 
of a single quasiparticle state to a change in the many-body distribution, and is thus 
non-zero only for interacting systems fH'. To leading order in any small expansion 
parameter, JF is equivalent to the matrix element for forward scattering ^T] . A review 
of the Fermi liquid picture applied to relativistic systems is given in jT^ . 

Departures from equilibrium are described by a transport (Boltzmann) equation 

dfl OOTl — * — * — * — * 

^ = ^ + V(5n.Vfceo-Vfcno.Vfe = /(n), (6) 

where we have used Hamilton's equations 

9^ ^ 9k ^ ,„, 

_ = V.. ; ^ = -Ve, (7) 

and dropped terms of 0{6^). The term I{n) on the right hand side is the collision 
integral describing the net number of scattering events into the phase space element 
d'^xd'^k per unit time. For oscillation frequency u) with mean free time r between 



quasiparticle collisions, the dimensionless parameter ujt governs which kind of exci- 
tation dominates. For ujt ^ 1 collisions reestablish local thermodynamic equilib- 
rium in each volume element, and ordinary first sound waves propagate with velocity 
(3i = ^Jd^plde ~ kp/ydjji. For ujt ^ 1, collisions are unimportant, and local thermo- 
dynamic equilibrium no longer holds. Since for a Fermi liquid r oc T~^, this situation 
holds near temperature zero, and the consequent excitations are known as zero sound. 
With T = J = 0, we have 

— * —* — * 

Vfe^o = (^prh ; V^no = -rh(3F5{e - /i) = -rh6{\k\ - kp), (8) 

where m.k = \k\, m.rh = 1. Specialising to oscillatory excitations at the Fermi surface 
of the form 6n = (5(£:— /i)$(m)e**^^'^~'^*\ we find the self-consistent zero sound equation 

i 

{uj - j3Frh.k)^{rri) = m.kSdkphr i dQJ='{rh, m')$(m')- (9) 

where the integral over the unit sphere is normalised to f dVt = 1, and the geometrical 
factor Sd = 2/(47r)ir(f ). 

3 Analytic Solution 

For the GN2+1 model considered here, the Fermi surface is a circle, and the Fermi 
liquid interaction a simple function of the angle 6 between quasiparticle momenta. 
Quasiparticle interactions in the GN model are mediated by a scalar boson a: to 
leading order in 1/A^/ in the chiral limit the interaction arises from the exchange 
contribution to forward scattering, is repulsive, and has been determined to be jH] 

^i(^) = TW7^ r(l-cos^), (10) 

where fic is the critical chemical potential at which the transition from chirally-broken 
vacuum to chirally-symmetric quark matter takes place at T = 0. To leading order 
in l/Nf fie = So, the dynamically-generated quark mass at T = /i = 0. For reasons 
to be made clear below, we will generalise (|1(J|) to read JF oc i? — cos 6 with R > 1 ^. 
Using (fTn|) in we derive 

/jnr 
—{R-cos{e-e'))<!>{e'), (ii) 



^In principle this could arise through a small bare quark mass tti, whereupon the exchange term 
yields R = 1 + 77z^//x^ p^. Note, though, that for to 7^ there is also an attractive contribution 
from the direct term, which moreover is no longer 0{1/Nf). 



The variable ( = s + ij is the ratio of phase velocity jSq = uj/\k\ to Fermi velocity 
Pf- The constant G parametrising the strength of the Fermi liquid interaction will 
be specified below. 

To find a solution we attempt an ansatz of the form 

^fn\ _Sr- ^" ^Q^" ^ + Cn sin 9 cos" 9. 
^ ' ~ 2^ C - cos ^ ^ ' 

n 

After substituting this into (fTTj) the integrals over 9' can be performed using 

d9 cosn9 2a" fd9sm9cosn9 

(13) 



2'7r(^ — COS0 a — a^ J 27r ( — cos 9 

where a, a~^ are the roots of z"^ — 2(z + 1=0 and we specify \a\ < 1. Equating 
coefficients of cos" 9 it is straightforward to show that the only non-trivial coefficients 
in (jl2|) are 6i and 62; we arrive at the coupled system 



hi = —[2610 + 62(1 + 0^ 

a — a '- 



62 = 



a — a ^ 



6i(l + a2) + ^(3a + a3) 



and immediately deduce 



69 a-2 - 1 - 2GR 



61 GR{a + a-i) 
Substituting this back into (fT^ determines a cubic equation for a^: 



(14) 



(15) 



,2 i\/r^„4 , „2/r,Dr^2 



l)(Ga^ + a'(2i?G' + G(3 - 4i?) - 2) + 2) = 0. (16) 

Discarding the trivial solution a^ = 1 we find a regular real solution with a < 1 exists 
for R > 2/(2 — G) and is given by 



1 



a^ = 1 + (2i? - |)G - G'R - V (1 + (2i? - |)G - Gmy -2G] . (17) 



^ \ ' ■ V— 2^ 

The corresponding sound speed is given by 

s = -(a + a-^)>l. (18) 

Note that since s > 1 the denominators of ()12ll3p are never zero; there is no 
necessity to specify a pole prescription to evaluate the integral and 7 = 0. Physically 
this means the wave is undamped. There are no solutions of the form (fT^ with 
complex a which describe dissipative solutions with s < 1 and 7 7^ 0. 
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We now return to the issue of why it is appropriate to consider R > 1 for the 
GN model in the chiral hniit. The expression ()1()|1 is derived from the relation 
J-" (6) = —A{6), where A is the forward scattering amplitude for physical quasi- 
particles on the Fermi surface, and is valid only in the large-A^j limit. Beyond 
this leading order it is necessary to consider a more general amplitude A4 {6, q) 
with momentum transfer q = {qo,q). For quasiparticle forward scattering at the 
Fermi surface A = lim|q^^oliKigo^o-^(9)) whereas for the Fermi liquid interaction 
JF = — limq^^^Qlim^^^Q Ai{q) j8j (see also JHl)- Taking the non-commutativity of the 
limits into account yields the relation 



Ai9) = HO) + 1^ / ^^(^')-^(^ - ^') 



(19) 



where g is the degeneracy of single particle states resulting from the trace in 0. Ex- 
panding A{6) = — J^n^n cos n9, J-'{9) = J2nQn cos n9, we readily derive the general 
relation 

Q" = ,_^p - (20) 

For the GN model the only non- vanishing coefficients are n = 0, 1; we recover ()11|) 
with 

G^^^G^^f';^ ,^,,,^^ ,+OiNj^y, R^'-±i> ' 



2 + G' SNfpFii^-i^,) 8iV/(/x-/Xe) ^ ^ ^' 2-G 2-G" 

(21) 

A zero sound solution with s > 1 therefore exists for all /i ^ IGNffic/i^QNj — g). 

In Fig. HI we plot s(/i) for the model defined hj q = 2, Nf = 4, ^cCt = 0.16 
corresponding to the simulation results to be presented in the next section. Note that 
chemical potential is expressed in cutoff units. The plot only shows those solutions 
with /ia ^ 0.18 for which the physical sound speed Pq = sjSp is definitely sub- luminal, 
using the relation lim^^oo /?f = 1 — b/32A^/ = || (note, though, that for finite /i this 
probably over-estimates j3f) |^- The plot shows s to be a rapidly decreasing function 
of /^, although the scale does not permit the observation that lim^_^oo(s(Ai) — 1) ~ 
0.202 X 10-^ > 0. 

To illustrate the form of the solution we plot ^{9) for various n in Fig. |21 The 
wavefunction has nodes at 9 = ±|, and is highly peaked in the forward direction, 
once again due to the causality restriction that s ~ 1 implying the denominator of 
(fT^ is very small as ^ — > 0. This would appear to be a general feature of weakly- 
interacting relativistic systems: indeed. Fig. |2l shows the peak sharpen as /i increases 
and hence G{fx) decreases, until a constant profile is reached as /x — >^ oo. 
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Figure 1: The ratio s = (3q/ (3f as a function of ^a for the GN model with Nf = 4 

4 Numerical Results 



We have simulated lattice-regularised versions of the GN model in 2+ld with both 
Z2 and SU(2)(8)SU(2) global chiral symmetries. The underlying lattice action uses 
staggered fermions: further details of the formulation and simulation algorithm are 
given in [21 E] • Most of our results are from the Z2 version, for which unless otherwise 
stated the simulation parameters are identical to those of our previous study [U], 
namely Nf = 4 and a/g"^ = 0.75, corresponding to a physical fermion mass at yU = 
of Eoa ~ 0.17. We have studied both 32^ x 48 lattices (the volume used exclusively 
in in]) and 48^; the increased momentum resolution offered by the latter has proved 
to be important. All results are taken in the chirally-restored phase, ie. with /i > 
/ic ~ 0.16a~^. Measured values for the baryon density in lattice units a = 1 are given 
in Table m and are to be compared to the saturation value ns(/x — ^ 00) = Nf/2. 

By analysing the decay of the appropriate correlation function with Euclidean 
time, we have calculated the dispersion relations E{k) for the spin-i quasiparticle, 
which carries a baryon charge, and uj{k) for various meson states of the form -^r-^ 
which are best thought of as excitations formed from a particle-hole pair. We have 
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Figure 2: Zero sound wavefunction $(^) for representative values of fia for the GN 
model with Nf = 4. 

simulated L^ x 48 systems with Lg = 32 and 48, with /i ranging from 0.2 to 0.6, and 
have measured E{k),uj{k) for k = {k,0) with k = 0,271/ Lg, . . . ,n/2. Note that the 
staggered fermion action is only invariant under translations of two lattice spacings, 
which restricts the space of accessible momenta. Further details can be found in ^. 
The quasiparticle propagator Q{k,t) is measured using the formalism outlined in 
|3], although in this case an external diquark source term is not introduced, so the 
"anomalous" propagator of the form {ip{0)ip^{x)) vanishes. In the chirally-restored 
phase Q{k,t) vanishes for t even; on the odd timeslices we fit to the form 



g{k,t) = Ae-^'^^''^' + Be~^^^^'^^^'~'\ 



(22) 



On 48^ acceptable fits were found for the most part for t G [7,43], although the 
fitting window had to be reduced somewhat for /i = 0.2, 0.3. Empirically we find that 
for small k the fit is dominated by a forward-propagating state (ie. A ^ B) which 
is interpreted as a hole, and for large A; by a backward-propagating state {A <^ B) 
interpreted as a particle. The transition between particle and hole behaviour is rather 
sudden, characteristic of a well-defined Fermi surface. In Fig. 121 we show E{k) with 
hole energies Eh plotted negative. Also shown are fits to the data with A; > of the 
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Table 1: Baryon density Ub = {'iplo'^') for the Z2 model with g ^ = 0.75. 



form 



E{k) = -Eo + Dsmh'^{smk), 



(23) 



ie. that of a Fermi liquid with discretisation effects taken into account, with effective 
Fermi momentum Kp and Fermi velocity Bp given by jl] 



Kp = sinh (sin kp) 



15 



Bf 



dsmh{E + Eo 



dsink 



D 



cosh Ei 







E=0 



cosh Kp ' 



(24) 



where E{kp) = 0. The resulting Kp, Bp and Kp/fiBp are tabulated in Table El and 
should be compared with the continuum values predicted in the large- A^/ approach 
for /i :» /ic 0: 



/3f 







0.984 



kf 



9 



0.969. 



(25) 



32Nf ' fipp 16 Nf 

The agreement is at best qualitative; the computed 0{1/Nf) corrections are so small 



/^ 


Kp 


Bp 


Kp/fiBp 


0.2 


0.190(1) 


0.989(1) 


0.962(5) 


0.3 


0.291(1) 


1.018(1) 


0.952(4) 


0.4 


0.389(1) 


0.999(1) 


0.973(1) 


0.5 


0.485(1) 


0.980(1) 


0.990(2) 


0.6 


0.584(1) 


0.973(1) 


1.001(2) 



Table 2: Fermi liquid parameters resulting from fits of (J23|l to data from a 48^ lattice 
(quoted errors are purely statistical). 

that much better control over systematic effects would be required for a meaningful 
comparison with the data. For instance, the Fermi liquid parameters should be most 
accurately pinned down for fia ^ 0.7 implying roughly equal numbers of points with 
k < kp and k > kp. Note too that analytic prediction of Pp, kp to 0{1/Nf) for finite 
/i strictly requires a knowledge of the gap equation to two loops, which is not yet 
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Figure 3: Quasiparticle dispersion relation E{}i) from simulations of the Z^ GN model 
on 48^. 

available [Hj. For the remaining analysis we will therefore assume the free field values 
Kp = /i, Bp = 1. 

We have observed a similar situation in simulations of the SU(2)®SU(2) model on 
32^ X 48 with coupling g'"^ = 0.85 chosen so that the vacuum fermion mass Sq matches 
that of the Z2 model at g^^ = 0.75. A much larger discrepancy in this model between 
the measured Fermi liquid parameters and the prediction (|25|) was reported in |1], 
probably due to the relatively small ratio {^ — ^1^/ ^i used in that study, implying 
that the Fermi liquid interaction is strong and higher order effects are important. 

Next we consider the meson sector, ie. correlators of the form 



Crik.t) = ^(^r^(0,0) ^r^(f,t))e^ 



ik.x 



(26) 



which carry zero baryon number. In Ref. [Bj we showed that the generic large-distance 
behaviour in any given channel is dominated by zero-energy particle-hole pairs, and 
as a result the decay is algebraic, ie. Cr{t) oc t~'^^^\ with the exponent A determined 
by the geometry of the overlap between two Fermi disks with relative displacement k 
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Figure 4: Meson propagators |C||| and |C_l| for both k 
The dashed hne shows a fit of the form (j28|) . 



and 



n 



f for /i = 0.5. 



between their centres. A particularly favourable configuration occurs for \k\ ^ 2/i, in 
which case the Fermi surfaces just kiss, and A = |. In this paper we focus on point- 
split meson operators corresponding to spatial components of the conserved vector 
current, ie. of the form 



Ji[X) 



r]i[x) 



[x(a;)x(a; + i)+ xi^ + i)xix)] , 



(27) 



where X:X ^^^ staggered fermion fields, ?7i(x) = (—1)*, and ri2{x) = (— 1)*"'''^^ Just 
as before we set k = (fc,0), and define C\\ in terms of the correlator (ji(O)ji(x)) 
and Cj_ ~ {J2{0)J2{x)). In Fig. HI we show |C||^j_|(/i;, t) data taken on a 48^ lattice at 
/i = 0.5, for both k = ^ ^ 2fi and fc = vr — | (the modulus is taken because we use a 
log-linear scale, and C fiuctuates in sign). Whereas Cj_ and C\\{k = |) all show the 
expected algebraic decay, the decay in the C||(fc = vr — |) channel is much faster, and 
resembles the exponential decay expected of an isolated pole. We have fitted it with 
the form 

C||(7r -k,t)= Aexp{-n{k)t) + (-l)*5exp(-cu(A;)t), (28) 
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in most cases employing datapoints with t G [10, 38]. 
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Figure 5: Dispersion relation uj{k) for various fi for the channel defined by the corre- 
lator C||(7r — k). Data are from 48^. 

For small k, the coefficient B ^ A, suggesting that the correlator is dominated by 
a pole in the alternating channel. In Fig. El we plot the resulting dispersion relation 
uj{k) for fi = 0.2, . . . , 0.6. The behaviour u oc k a.s k ^ suggests the presence of a 
massless pole similar to that of a phonon. Support for this interpretation comes from 
recasting the meson bilinear in terms of continuum- like fields g"°, g"° having spinor 



index a = 1, ... 4 and 'color' index a = 1, 2 



We obtain 



1)"^ {-iyx{x)viX{x + 1) ~ tqilo ® r*)g. 



(29) 



demonstrating that the excitation can be viewed as an oscillation of local baryon 
density. 

Since we are attempting to simulate low temperatures T <^ /i, our expectation is 
that the phonons are characteristic of zero sound and should be regarded as collective 
excitations of the degenerate ground state, rather than of first sound which relies on 
particle collisions to generate the elasticity of the medium, and which is therefore 
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heavily damped for small T ^Uj. An obvious objection to this proposal is that the 
phase velocity Pq = u/k derived from Fig. El has Po ~ 0.55 - 0.65, whereas the 
theoretical considerations of Sec. El imply that undamped oscillations in relativistic 
matter should have /5o ~ 1. However, recall that although we are examining the 
k ^ limit, we are trying to model a Fermi surface phenomenon on the lattice; the 
appropriate scale with which to compare /3o is the "bare" Fermi velocity (ie. without 
the discretisation correction modelled by (|23I24|) . which for free fields is given by 



/obare 

Pf 



BE 
'dk 



1 — sinh^ jj, 



E=o \ 1 + sinh /x ' 



(30) 
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Figure 6: Sound velocity /3o rescaled by the bare Fermi velocity Pp^'^'^ for various /i. 

In Fig. iniwe plot Po/ 0^^^ for fc > using the 48^ data of Fig. O There are two 
striking features. First, the ratio appears to extrapolate to a value rather closer to, 
though still smaller than, 1 as A; ^ 0, getting closer as /i increases until for yU = 0.6 
the intercept > 0.95. Secondly, for A; < || the curves are remarkably close to straight 
lines, implying a dispersion of the form uj = sk — Fk"^. Fig. IHlalso shows fits of this 
form, and the resulting intercept and slope are given for all the systems studied in 
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TableEl The most successful fits were to data on 48^ with /i = 0.6, where straight hue 



1^ 


48^ 

s F 


322 X 4g 
s F 


Z2 {g-^ = 0.75) 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 


0.896(8) 0.578(4) 
0.903(5) 0.576(20) 
0.911(4) 0.500(16) 
0.924(3) 0.447(9) 
0.959(2) 0.472(4) 


0.894(5) 0.552(15) 
0.961(3) 0.575(8) 
0.950(2) 0.425(2) 
1.033(3) 0.556(8) 


Z2 ig-' = 0.6) 
0.6 


0.977(2) 0.460(5) 


_ _ 


SU(2)®SU(2) {g-^ = 0.85) 
0.5 Is = 1 

0.5 Jg = 


- 


0.938(2) 0.467(7) 
0.917(3) 0.453(9) 



Table 3: Straight line fits to uj/kj3p^ vs. k (quoted errors are purely statistical). 

fits to data with A; < ^ proved possible. For /i < 0.4 the fits only included data with 
fc < I (48^) or fc < fl (32^ x 48). The general trend revealed in Table El is that the 
speed ratio s increases towards unity as /x increases, and that the slope parameter F 
decreases slightly. Since the momentum resolution on 32^ x 48 is markedly worse, we 
can only comment on the absence of evidence for any significant finite volume effects, 
and, in the model with global SU(2)®SU(2) symmetry, for any significant difference 
between isosinglet ('u;') and isotriplet ('p') channels. 

The most interesting comparison is between Z2 data on 48^ with g'"^ = 0.6 and 
0.75, the stronger coupling corresponding to roughly a factor of two increase in lattice 
spacing; if we restore explicit factors of a we see that the slope parameter Fa~^ 
remains roughly constant under this change, implying that the slope is a lattice 
artifact which should vanish in the continuum limit. At first sight, it is a little 
strange to observe an 0{a) scaling violation in a simulation with staggered fermions. 
However, a careful analysis |T3] reveals that there are indeed 0{a) discretisation 
artifacts in the lattice formulation of the GN interaction term {^ipY. This highlights 
that the effect we are examining is a property of an interacting system, since non- 
interacting staggered fermions show discretisation effects only at 0{a?). Since s also 
increases as the coupling is strengthened we deduce a systematic effect as ^a and 
hence p/T increase. 
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5 Discussion 

We have, in Sec. El used the Fermi hquid interaction to leading non-trivial order in 
1/A^/ to find solutions to the Boltzmann equation corresponding to zero sound. The 
excitation is spin and isospin symmetric, and has speed s = (3q/ (3p > 1 for almost 
all /i > fie- In the simulations described in Sec. El we have found a massless pole in 
a particular mesonic channel, interpolated by a staggered fermion bilinear point-split 
along the direction of propagation, with continuum quantum numbers corresponding 
to a local fiuctuation in baryon density, with speed s ^ 1. To what extent can we be 
sure these results describe the same physical phenomenon? 

Our solution based on the ansatz (fT^ . is the unique one (symmetric in internal 
quantum numbers) which is a meromorphic function oi z = e*^. As remarked above, 
the poles fall on the real axis, yielding a real s > 1. Solutions with s < 1 must 
therefore be of a different form to ()12j) . and probably involve a branch cut ending 
inside the |2;| = 1 contour. We can see this physically by considering the emission of 
a phonon with momentum g'and energy l3o\q\ from a quasiparticle with momentum 
k and energy fi + Pp{\k\ — kp). This is allowed kinematically for s < 1, the angle of 
emission (p satisfying 

COS(l) = S+Mr(l-S^)> S. (31) 

All radiation is emitted within a cone of half-angle cos~^ s centred on the quasiparticle 
trajectory. This well-known phenomenon is variously known as Landau damping, 
Cerenkov radiation, or most appropriately in the current context, as a sonic boom. 

On a spacetime lattice, this process is constrained because there is a natural 
lower bound for the angle of emission, ^min ~ 27r/Ls|g|. Landau damping is thus 
kinematically forbidden for 

S > COS0min ^ 1 - T7T-q^- (32) 

With a conservatively high \q\ ~ fi/2, then on a 48^ lattice we have no damping 
for s > 0.73 at /i = 0.2, rising to s > 0.97 at /i = 0.6. It is plausible therefore that 
damping is suppressed and that the phonon found in Sec. El is described by an isolated 
pole even if s < 1; indeed. Fig. IHl suggests that s increases towards one as fi increases. 
This contradicts the analytic solution shown in Fig. Q which has lim^_^oo s -^ 1+; 
note, however, the disparity in vertical scale between the two figures. Simulations on 
volumes considerably larger than those used here will be needed to disentangle the 
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various possible systematic effects due to finite L^, finite /x/T, and non-zero a in order 
to determine the sign of s — 1 for finite Nf. In principle, since our extracted value of 
s ^ 1 has relied on a rescaling of jSp to take account of discretisation artifacts, we also 
need to go to considerably finer lattices in order to demonstrate a clear distinction 
between our signal and first sound with expected propagation speed /?i ~ 0.7. If, 
however, the theoretical arguments about zero sound dominating as T ^ can be 
taken seriously, then for the first time we have succeeded in identifying a collective 
oscillation in a lattice simulation. 
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